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COST PERFORMANCE SATELLITE DESIGN USING QUEU NG THEORY 

Gerald F. Hein 

Space Flight Programs Directorate 
NASA Lewis Research Center 
Cleveland, Ohio 


ABSTRACT 

The use of satellites for wide area direct broadcasts results in 
large power requirements for the satellite. Such satellites are or will 
be used to provide remote health care delivery nnd educational services 
to remote regions and many other new services including warnings to the 
public of impending natural disasters. These new special purpose satel- 
lites are very cost sensitive to the number of broadcast channels, usually 
will have Poisson arrivals, fairly low utilization (less than 35%), and a 
very high availability requirement. To solve the problem of determining 
the effects of limiting C the number of channels, the Poisson arrival, 
infinite server queu'ng model will be modified to describe the many 
server case. The model is predicated on the reproductive property of the 
Poisson distribution. 

For smali changes in the Poisson parameter under the assumptions 
stated previously, the resulting distribution of states or number in the 
system will be Poisson. A difference equation has been developed to de- 
scribe the change in the Poisson parameter. When all initially delayed 
arrivals reenter the system, a (C + 1) order polynomial must be solved to 
determine the new or effective value of the Poisson parameter. When less 
than 100 percent of the arrivals reenter the system, the effective value 
must be determined by solving a transcendental equation. 

The model was used to determine the effects of limiting the number 
of channels for a disaster warning satellite. State probabilities and 
delay probabilities were estimated for several values of the number of 
channels C for arrival and service rates obtained from disaster warn- 
ings issued by the National Weather Service. The results predicted by 
the queuing model were compared with the results of a digital computer 
simulation. 


INTRODUCTION 

The use of geosynchronous orbit satellites for communications began 
in the 1960's with SYNCOM I. Since then a network has been established 
to provide international telecommunications vja satellite to about 80 
countries throughout the world. The use of modem telecommunications has 
become so commonplace that there is often a tendency to forget how recent 
are the technological achievements which make such usage possible. 

The present network for commercial traffic is used for point -to- 
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point transmission of telephone and television. Such a network requires 
large and expensive ground transmitter-receiver systems. 

In late 1971 an initiative study was conducted for the Executive 
Office of the President primarily by NASA and included about 100 partici- 
pants from various government agencies. This study, "Communications for 
Social Needs," was an effort to solve social problems through the use of 
advanced technology. One of the results of the study was the conclusion 
that there is an evolving need for what might be called special purpose 
communications satellites. 

The initiative study concentrated on an area of communication satel- 
lite applications In which the ground receiver systems or terminals were 
very large in number and hence a driving cost parameter. Contrasting the 
commercial and the special purpose types of communication satellites , one 
finds the following differences: 

Commercial : Low power, point-to-point transmission, small number of 

expensive terminals, large number of channels 

Special purpose : High power, wide area coverage, large number of 

low cost terminals, small number of channels 

In usage, the operation of the special purpose type generally con- 
sists of transmitting information to many receivers over a relatively 
large geographic area. The applications considered in the aforementioned 
study were the uses of communications for remote health care delivery, 
electronic mail handling, law enforcement, education, and as a possible 
means of warning the public of Impending natural disasters, such as hurri- 
canes or tornados. 


The most important difference between the commercial and the special 
purpose satellites is that the major design objectives are radically op- 
posed to one another. Herbert Raymond^ suggested that the most meaningful 
parameter for a commercial satellite system is probably the utilization 
factor. It should be high In order to maintain a profitable system. The 
special purpose satellite must be designed using availability as the major 
parameter. In the remote regions of the far northern hemisphere, for ex- 
ample, where neighboring villages may be separated by distances of hun- 
dreds of miles , the need for emergency medical care can be met through 
the use of paraprofessionals communicating with doctors via satellite. 
However, it is not reasonable to expect anyone to wait for service in 
this application. Should satellites be used in the future to warn the 
public of impending natural disasters, it would be essential to dispatch 
warnings as quickly as possible. A system which warns of tornados 15 min- 
utes after their occurrence would be of little value to the public. 


The National Aeronautics and Space Administration (NASA) and the 
National Oceanic and Atmospheric Administration (NOAA) have been conduct- 
ing joint investigations of various technologies in order to examine the 
feasibility of using communication satellites for one of the applications 
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mentioned previously, namely, to provide warnings to the general public in 
the event of an impending natural disaster. The various candidate systems 
for disaster warning which have been suggested for consideration include 
the mass ringing of telephones, microwave transmission of radio signals, 
terrestrial radio networks and the use of communication satellites. 

Government organizations other than NASA are conducting studies of 
terrestrial systems and NASA is confining its investigation to the use of 
satellites. When completed, the studies will be used to determine the 
most cost effective system. At the present time, satellites offer a very 
viable alternative because several meteorological functions may be com- 
bined with the communication function, and the satellite system has the 
desirable property of being "hardened" against natural disasters. That 
is , satellites are not prone to destruction from an impending natural 
disaster. 

The functions of a natural disaster warning system as reported by 
Hein and Stevenson^ are: 

1. Route disaster warnings to the general public 

2. Provide disaster communications among national, regional, and 

local weather offices and affected areas 

3. Provide environmental information to the public 

4. Provide a system for collecting decision information for the dis- 

semination of warnings 

The natural disasters which would be monitored by a DWS include tornados, 
severe thunderstorms, flash floods, tsunami, earthquakes, hurricanes, 
forest fires, winter storms, and a category called other. 

In September 1973, the Aerospace Systems Group of the Computer 
Sciences Corporation in Arlington, Virginia, began work on a feasibility 
study of using satellites for a natural disaster warning system. This 
study, funded by the NOAA and managed by the NASA Lewis Research Center, 
provided the conceptual design of a satellite DWS. 

As will be indicated in more detail, the primary problem with the de- 
sign of such a satellite DWS is the determination of the minimum number of 
communication channels required to service the system needs. This deter- 
mination requires an extensive analysis of the message traffic for the 
proposed system in the early 1980's. 

The satellite system will require a high Effective Isotropic Radi- 
ated Power (EIRP) . Present satellite power is on the order of 50 watts 
per channel and designs of the mid- to late seventies will be limited to 
less than 500 watts per channel. As reported by the Computer Sciences 
Corporation, DWS satellites will require more power than present designs 
can deliver. If the satellite is to ba a viable alternative, the total 
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power will be limited to something less than 10 kilowatts which thus pre- 
cludes the addition of large safety margins for the channel requirements. 

There is a genuine need to determine the minimum number of channels 
for the DWS requirements. The application of satellite technology for 
the solution of social problems such as those discussed in the study "Com- 
munications for Social Needs" may be delayed if the problem discussed 
above is not solved. The requirements for the DWS satellite are very sim- 
ilar to the requirements for a large class of satellites which may be 
used to solve some social problems. For lack of a better description, 
this class of satellites will be called special purpose communication 
satellites. The salient need is to determine the performance criteria 
required for the estimation of the channel requirements. The alternatives 
are to overdesign which is not feasible here, use simulation techniques 
which may be expensive or use an analytical technique. 

The potential solution of the channel requirements problem for spe- 
cial purpose communication satellites was thought to be in the problem 
class of the GI/G/C queue. An extensive analysis was made of the types 
of messages sent over the teletype network of the National Weather Ser- 
vice. It was found that the messages could be classified into six differ- 
ent groups of input. It was also demonstrated that the six types of mes- 
sages could each be represented by a Poisson distribution. Although such 
a simplification should be expected for large numbers of messages over 
the period of a month, the hypothesis testing was done for relatively 
short intervals such as four or five days because of requests by those 
involved with the problem at NOAA. It is very important to know the dis- 
tributions of what NOAA personnel have called "spikes" or short bursts 
of increased traffic Intensity such as that experienced during the rampage 
by Hurricane Agnes from June 14 to 28, 1972 or the tornados of April 3, 
1974, so that the communication system can be designed to handle such 
traffic loads. 


ANALYTICAL MODEL 

To develop the framework of the model, reference will be made to the 
Poisson arrival model with an infinite number of servers , The Infinite 
server case has been solved for Poisson arrivals. The basic hypothesis 
is that the state distribution for the many server case should be similar 
to that for the infinite server case. Because of the reproductive prop- 
erty of the Poisson distribution and the fact that the mean and variance 
are equal, a change in the Poisson parameter X to X + AX results in a 
Poisson distribution with a mean and variance of X + AX. 

In the case of an infinite number of servers with Poisson arrivals, 
Takacs proved that with an arbitrary service discipline, the limiting dis- 
tribution of the states can be represented by a Poisson distribution. 
Moreover, the distribution is independent of the Initial state. The 
proof of this theorem is given in Takacs^ and is stated here. 
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Let 



X dH(X) and 


P[5(t) = K] 


P^Ct) where a is the average 


service time, P{ t (t) is the probability of being in state K at time t, 
0 is the Poisson arrival parameter, and £(t) is the queue size or num- 
ber in the system. 


Theorem 1 . If 5(0) = 0, then 



( 1 ) 


lim P K (t) = P* (K-0, 1, 2, ...) 

fc-Ko ^ 


exists and we have 


P * . .-*■ 

r K k: 


( 2 ) 


In the case of the many server queue without storage, arrivals will have 
a certain probability of being rejected if there is no server available 
when the arrival enters the system. In the case of the special purpose 
communication satellite, most of the arrivals will keep trying to enter 
the system until a server is available. In many problems, however, a 
fixed portion will leave the system and not return. The fixed fraction 
not returning will be specified by (1 - y) where 0 <_ y <_ 1, 

As stated previously, y will usually be 1 in the application here. 

In the theorem of Takacs , the state probabilities of the system can 
be described with a poisson distribution. The parameter for this distri- 
bution is the Poisson parameter for the arrival rate times the average 
service time. That is , A = 8a. When the number of servers is con- 
strained, the effect of the constraint may be described as an effective 
increase in the average number of arrivals. Since the service time re- 
mains constant, a difference equation can be written to describe the ef- 
fects of server -customer interaction. An arbitrary time Interval will be 
selected so that the rejected arrivals will have a very low probability of 
entering the system more than once during the same time interval. For the 
first interval of time at startup, the parameter X± - d±a may be de- 
scribed as 


1 •. 


X 0 = X 0 + 0 


During the second time interval 


X x - X Q + YX 0 4>(C,X 0 ) 


where y is the portion of rejects returning; $ is the complementary 
cumulative of the Poisson distribution; Xq is the Poisson parameter, and 
C is the number of servers in the system. 

To simplify the notation, (i) will be used to denote 4>(C,Xi). In 
general, 

X. - X Q [1 + Y4>(.1 - 1) + Y 2 $(j - l)*(j - 2) + . . . + Y J <Kj - 1) . . .*(0)] 


= x Q ti + 


lid: H-l 1 


The proof of convergence for this sequence is given in Hein.' 
C = n, Y ° 1 . For the general n server case 


“ x * , . .2 , , * . * 

X A = e I X , + X , + . . .+ 7 t-TT + — i — 

0 r * * (n - 1 ) . n! J 

Differentiating with respect to X A and equating to zero yields 


"(n - 1) ! 


• - ** - 1 - 0 


According to the rule of signs developed by Descartes, this polynomial 
has exactly one positive root. To obtain the roots of polynomials there 
are many computational methods available. One which is rather easy to 
use is the Newton-Raphson method. According to this method the next 
estimate of the root of f(x) is given by 

f(x n ) 

X n+1 X n f ' (x ) 
n 


After the value of is obtained Xqj,j^ may be computed. If 

XoMAX i s less than C, then Xgj^y; is the maximum ),q for that C; 
otherwise C is the limit. 
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APPLICATION OF MODEL 

An extensive analysis of six years of NWS message traffic demon- 
strated that message arrivals in the five categories : hurricanes , tor- 

nados , winter storms, small craft, river and other could be described 
using Poisson distributions. Broadcasting times for all categories ex- 
cept hurricanes may be described using uniform densities. Hurricane 
message broadcasts may be described using a log normal density. Complete 
results of the statistical analysis are reported in Hein.l 

Using time units of one minute, estimates of relative message fre- 
quencies and broadcast times for December were calculated from NWS data, 
December was selected because it yields the highest traffic density and 
the longest average broadcast times. 

The Poisson parameter and the average broadcasting time may be used 
to determine the initial state parameter: 

X 0 = Q Q a 

Once the value of Xq has been determined, it is necessary to determine 
the feasible range for the number of channels. If Xq is less than 
XoMAX for one channel, the range of feasibility is 1 <_ C <_ <=°. The 
equation 



must be solved for X^ over the range of values of interest for C, the 
number of simultaneous broadcast channels in the satellite. After deter- 
mining X A for the range of C, the state probabilities must be deter- 
mined from the equation 


k 


p 


K 



K = 0, 1, 2, . . . 


These state probabilities for the effective lambda (X^.) may be used to 
calculate the probability of a delay, or to determine the probability 
of a delay exceeding a given time value. Now the determination of the 
effects of limiting C becomes a relatively simple analytical procedure 
and does not require the use of simulation. The procedure may be used to 
determine the required number of channels for any queuing system with 
Poisson arrivals and arbitrary service. A criterion may be established 
in advance and then the implications of the criterion may be analyzed. 

In order to determine the delay or queuing time for messages, it is 
necessary to introduce a time factor, The probability of a delay is the 
probability that Wq (queuing time) is greater than zero: 
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P(W > 0) 

q 



P 

n 


where C Is the number of simultaneously accessible channels in the 
satellite system. The average service rate per channel is 1/a. For 
C channels this rate is C/a. But this is true only when there are C 
channels available fo: service. 


The average service rate is then 



The complementary waiting time distribution is then given by 


P{W q 


t} 


e -C/a(l-X*/C)t p £ W > 

q 


SIMULATION RESULTS 

The basic time unit for the process reported here is one minute. Be- 
cause the basic time unit is small, many problems normally associated with 
obtaining a steady state have been avoided. The parameter of prime im- 
portance here is the number of message requests in the system at any given 
time. Short simulation runs of one or two days were more than adequate 
to demonstrate that the process cycles "Infinitely often" through a single 
state (namely zero). When such a cycling occurs, the sample sizes re- 
quired to reduce intervals of uncertainty are significantly smaller than 
what would normally be required because steady state is achieved upc \ 
startup. Extensive digital simulations were conducted to verify the ana- 
lytic model. The simulations and the statistical analyses of the results 
are thoroughly documented in Hein.-*- Copies of the computer software and 
sample results are also included. 


SUMMARY AND CONCLUSIONS 

The results obtained from the model provioed excellent approxima- 
tions for the 6- and 4-channel cases. Statistically there was no signif- 
icant difference between the predicted results and those obtained from 
simulations. When the traffic intensity was increased through a reduc- 
tion in service channels to 2, the customer-server interaction caused a 
degradation in the quality of the predicted results. The interaction 
causes the variance of the Poisson distribution to increase through a 
distortion of the relative state frequencies. This distortion causes a 
change in the state distribution; it deviates from the predicted Poisson 
distribution. Although the 2-channel case may provide useful results. It 
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also demonstrates that the iiodel is close to the limits of usefulness, 

As a result of the model described here, a more realistic appraisal 
of the channel requirements for a satellite DWS was obtained. There are 
many areas of application beyond usage of the model for communications 
satellites. Applications to queuing systems where low utilization/high 
availability is the predominant characteristic are more numerous than can 
be mentioned here. Moreover, the necessity of minimizing waiting time 
will become even more important in the future as world economies become 
dominated by service industries. Market strategies will include trade- 
off analyses of the value of minimizing waiting versus the cost of adding 
more service channels. 

The primary value of the model is that it 'educes the dependence on 
expensive computer simulations when availability is a prime criterion of 
design. Once verified, the model allows rapid calculation of system 
statistics and provides reliable information at relatively low cost. 
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ALTERNATIVES 

1. USE SIMULATION TECHNIQUES 

2. USE CURRENTLY AVAILABLE 
ANALYTICAL TECHNIQUES 


3. EXTEND CURRENT ANALYTICAL 
TECHNIQUES 
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Let a- / X dH(X) and P [C (t) - K] <= P k (t) 
'0 


where a is the average service time, P k (t) is c ^ e probability of 
being in state K at time t, 0 is the Poisson arrival parameter, 
and £(t) is the queue size or number in the system. 


Theorem 1. If £(0) = 0, then 


■ 


and if a < “, 




for K = 0, 1, 2, 


lim P k (t) = P^ 

t-+03 


(K = 0, 1, 2, . . .) 


exists and we have 


* _ -0a (_0a) K 

K k: 


x i - V 


x o = x o + 0 


x ! “ x 0 + YX 0 $(C,X 0 ) 


In general, 

Xj = X Q [1 + y$(j - 1) + y 2 $(j - 1) 4>(j - 2) + . . . + Y j *(j - 1) 


j=l 


•% 

V"* 

1-1 


1 + ) Y** 

"|~J *(K) 

> 

l—J 

_K“i _ 




J 


• • .*( 0 )] 
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r .-**/ X l X 

1 -e ll + X^ + 75-r + • • ,+ Q' 


* X l -h 

P K = K T e 


K = 0, 1, 2, . . 


n 


P(W q > 0) 



n=C+l 


P(w q > 0) - 1 - ; P n 

n=0 


n 


P(W q > t) = 2j? n PC w q > t|E n ) 
n=0 


P(W > t) = e 

q 


-C/a(l~X /C) t 


P (W > 0) 


and 


P(W <t) = 1 - P(W > t) 
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RATE TO 23,000 PER MONTH. 
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SUMMARY OF RESULTS 

ANALYTICAL TOOL IS RELIABLE AND SIMPLE TO USE. 

COST OF OBTAINING INFORMATION SIGNIFICANTLY REDUCED. 
COST OF SATELLITE SYSTEM SIGNIFICANTLY REDUCED. 


